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We consider a model where sterile neutrinos can propagate in a large compactified extra dimension 
giving rise to Kaluza-Klein (KK) modes and the standard model left-handed neutrinos are confined 
to a 4-dimensional spacetime brane. The KK modes mix with the standard neutrinos modifying 
their oscillation pattern. We examine former and current experiments such as CHOOZ, KamLAND, 
and MINOS to estimate the impact of the possible presence of such KK modes on the determination 
of the neutrino oscillation parameters and simultaneously obtain limits on the size of the largest 
extra dimension. We found that the presence of the KK modes does not essentially improve the 
quality of the fit compared to the case of the standard oscillation. By combining the results from 
CHOOZ, KamLAND and MINOS, in the limit of a vanishing lightest neutrino mass, we obtain 
the stronger bound on the size of the extra dimension as ~ 1.0(0.6) /im at 99% C. L. for normal 
(inverted) mass hierarchy. If the lightest neutrino mass turn out to be larger, 0.2 eV, for example, 
we obtain the bound ~ 0.1 /im. We also discuss the expected sensitivities on the size of the extra 
dimension for future experiments such as Double CHOOZ, T2K and NOi^A. 

PACS numbers: 14.60.Pq,14.60.St,13.15.+g 



I. INTRODUCTION 

Our observable 1+3-dimensional universe could be a 
surface, the brane, embedded in a dimensionally richer 
1+3+d-dimensional spacetime (d being the number of 
extra dimension), the bulk. This intriguing idea can be 
motivated by string theory, where at least 6 extra spa- 
tial dimensions are required, in particular, by stringy in- 
spired models designed to address the disparity between 
the electroweak (~ 1 TeV) and the gravity (~ 10 16 TeV) 
scales. There are two basic scenarios commonly evoked 
to generate the hierarchy between these two fundamental 
scales of nature: either by suggesting the source of the 
hierarchy to be the volume of a fiat extra dimensional 
space [T] or the strong curvature of that space 0. 

In this paper we are interested in constraining the large 
extra dimension (LED) scenario [1J in connection with 
neutrino physics since right handed neutrinos (standard 
model (SM) singlet fields) in this case can, as well as 
gravity, propagate in the bulk. Tablctop experiments de- 
vised to test for deviations of Newtonian gravity can only 
probe LED up to submillimeter sizes. The most stringent 
upper limit given by a torsion pendulum instrument is 
200 /xm at 95% C. L. for the size of the largest flat extra 
dimension regardless of the number of d [3]. Neutrino 
physics can be considerably more sensitive to LED. 
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We should, however, mention that astrophysical 
bounds on LED are in general much more stringent (see 
e.g., [I], and references therein) than the ones obtained 
by the terrestrial experiments including that from col- 
lides. However, these astrophysical bounds are not com- 
pletely model independent and therefore, we believe that 
studying the possible impact of LED which can be probed 
(independently from astrophysical constraints) by terres- 
trial experiments is still worthwhile. 

There are mounting evidences from several solar [5] , at- 
mospheric [5] and terrestrial [THT2"] neutrino experiments 
that neutrinos undergo flavor oscillations due to mass 
and mixing. As it was shown in [T51U7| . LED can have 
strong impact on neutrino oscillation probabilities. How- 
ever, since the current neutrino data mentioned above 
are perfectly consistent with the standard three flavor 
oscillation scheme, the effect of LED, if it exists, is ex- 
pected to be present only as a subdominant effect on top 
of the usual oscillation. Therefore, as was done in |18j . 
in this work we assume that LED effect would only per- 
turb somewhat the standard oscillation pattern and try 
to constrain LED using the current oscillation data. 

In this paper, we studied the possible impact of 
LED on the former and current oscillation experiments 
CHOOZ HU, KamLAND [HIE] and MINOS [MIO] in 
order to obtain the upper bound on the size of the largest 
extra dimension, which turns out to be submicrometer 
range. We do not consider solar and atmospheric neu- 
trino data in this work because the analysis would be- 
come much more complicated due to the matter effect 
and also because we expect similar bounds from these 
data (see Sec. M). We also calculate the expected sensi- 
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tivities on LED for future experiments such as Double 
CHOOZ [20], T2K [21], and NOz/A [221(23]. 

This paper is organized as follows. In Sec. [IT] we de- 
scribe the framework of our study of neutrino oscillations 
with LED. In Sec. IIIII and IIVI we discuss the best cur- 
rent and future limits that can be established on the size 
of the largest extra dimension from neutrino oscillation 
data. Finally Sec. [V] is devoted to discussions and gen- 
eral conclusions. In Appendix |A] we describe the solution 
of the neutrino evolution equation for a constant matter 
potential whereas in Appendix [B] we describe the details 
of our x 2 analysis. 

II. NEUTRINO OSCILLATION FORMALISM 
WITH LED 

We consider here the model discussed in Refs. [T6HT8] 
where the 3 standard model active left-handed neutrinos 
fields v^l (a — e, [i, r), as well as all the other SM fields, 
including the Higgs, are confined to propagate in a 4- 
dimensional brane, while 3 families of SM singlet fcrmion 



fields can propagate in a higher dimensional bulk, with at 
least two compactified extra dimensions (d> 2). We will 
assume that one of these extra dimensions is compactified 
on a circle of radius a, much larger than the size of the 
others so that we can in practice use a 5-dimensional 
treatment. 

By this assumption, our bounds are always more con- 
servative than the ones obtained by assuming all the LED 
radius a are the same for d > 2. In other words, if we 
have adopted the same assumption (of equal raduis for 
all LED), we should have obtained stronger bounds on a 
for d > 2 because the conversion into KK modes would 
be more efficient under such an assumption for d > 2. 

The 3 bulk fermions will have Yukawa couplings with 
the SM Higgs and the brane neutrinos ultimately leading 
to Dirac masses and mixings among active species and 
sterile Kaluza-Klein (KK) modes. The 4-dimensional La- 
grangian which describes the charged current (CC) inter- 
action of the brane neutrinos with the W as well as the 
mass term resulting from these couplings with the bulk 
fermions in the brane, after electroweak symmetry break- 
ing and dimensional reduction, can be written as [T8] 
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where the Greek indices a,/3 = e, fi, r, the capital Roman evolution Eq. (All I in the Appendix, and U and W are 
index N — 1, 2, 3, oo, m®p is a Dirac mass matrix, 

^oJi> u aR ano - v ctL are the linear combinations of the 
bulk fermion fields that couple to the SM neutrinos i/w . 

After performing unitary transformations in order to 
diagonalize m®o we arrive at the neutrino evolution equa- 



the mixing matrices for active and KK neutrino modes, 
respectively. 

This transition probabilities, even in vacuum, depend 
on the neutrino mass hierarchy since both W^ N) and 
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are functions of the dimensionless parameter 
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tion (A7) that can be solved to obtain the eigenvalues 
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(see Appendix |a| , so that 
the transition probability of into (subscript L is 
omitted) at a distance L from production, 



can be given in terms of the transition amplitude 
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where E is the neutrino energy, L is the baseline dis- 
tance, A!; is the eigenvalues of the Hamiltonian of the 



y/2m,j a, where mj (j = 1, 2, 3) are the neutrino masses. 
We will consider here two possibilities for the mass hier- 
archy: normal hierarchy (NH) with my, > rri2 > mi = tuq 
and inverted hierarchy (IH) with > mi > 772 3 = ttiq. 
As mo increases NH and IH become degenerate. We de- 
fine the mass squared differences as Amfj = m\ — mj 
(i,3= 1,2,3). 

To understand qualitatively the results to be presented 
in Sees. |III| and |IV| we discuss here what is to be expected 
of the effects of LED on the survival probabilities. In 
Fig. [T] we show the survival probabilities for and v e 
as a function of the neutrino energy E in vacuum for 
NH and IH for MINOS (735 km), KamLAND (180 km) 
and CHOOZ (1 km), assuming m = and a = 0.5 /im. 
There are three basic effects of LED: a displacement of 
the minima with respect to the standard survival prob- 
abilities, a global reduction of the flavor survival prob- 
abilities as SM neutrinos can oscillate into KK modes 
and the appearance of extra wiggles on the probability 
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FIG. 1: In the top (middle and bottom) panel we show the 
survival probability for f M (i/ e ) as a function of the neutrino 
energy E for the baseline L = 735 km (180 km and 1 km) 
for a = (no LED, black curve) and a = 0.5 fira for NH 
(blue dashed curve) and IH (red dotted curve). The other 
oscillation parameters were set to sin 2 612 = 0.32, sin 2 623 = 
0.5, sin 2 26>i3 = 0.07, Amli = 7.59 x 10~ 5 eV 2 , and |Amli| = 
2.46 x 10 -3 eV 2 . The lightest neutrino mass, mo, was set to 
zero. 



pattern due to the fast oscillations to these new massive 
modes. 

When matter effects can be ignored, the impact of LED 
in the survival amplitude, to leading o rder i n is such 

.(LED) 



that A y 



«Ei£ Itfaif (see Eq. 



A21). Therefore, 



roughly speaking, in order to modify the standard prob- 



ability by say ~ 10% by the effect due to LED, at least 
one of the £ 2 should be order of ~ 0.1. Since we know, 
from atmospheric neutrino oscillation, that at least one 
neutrino has a mass rrii ~ 0.05 eV, we can estimate that 
~ 0.1 implies a ~ 5cV _1 = 1 /im. So one can ex- 
pect the terrestrial experiments to be sensitive around 
this scale, which is consistent with our results discussed 
in the next section. 

From now on we will drop the (0) superscript when 
referring to flavor oscillations. In the case of — > v^, 
from Fig. [T] we see that the effect of LED is basically 
the same for NH and IH. This is because ^4 , D , f0) oc 

l^«| 2 is °f the same order for NH (mainly driven 
by £3) and IH (mainly driven by £2)- O n the other hand, 
in the case of v e — > v e the effect of LED is significantly 
larger for IH than NH since A_^° _ (0) oc \ U e t\ 2 = 

£i=i 2 £i I C^ei 1 2 + £3 sin 2 0i 3 is suppressed due to small 
sin 2 9x3 f° r NH (since £3 3> £1,^2 for vanishing mo) 
whereas for IH the dominant LED contributions (due to 
£1 and £2) are not suppressed. 

In this paper we do not consider the appearance chan- 
nels such as — > u e and — > P e due to the following 
reasons. First of all, when #13 is zero or much smaller 
than the current bound, we found that the impact of 
LED for these appearance channels is very small com- 
pared with that of the disappearance modes considered 
in this work. In principle, even if 0\3 is zero, LED can 
induce a flavor transition such as Vn — > v e (for T2K 
and NOi/A) through the right handed KK modes but 
such a transition is a kind of second order effect (this 
is because, ignoring oscillation driven by solar parame- 
ters, would not be converted directly to v e but only 
through KK modes) whereas the impact of LED for the 
disappearance channel is the first order effect, or it is 
the consequence of the direct transition from active to 
sterile KK modes. This argument applies also to the ap- 
pearance experiments like LSND [53], Karmen [25] and 
MiniBOONE and therefore we do not consider these 
experiments in this work, as they do not make any signif- 
icant contribution to improve the bounds we obtained. 

On the other hand, if &\3 is large enough to be ob- 
served by T2K and NOz/A (in the absence of LED), then 
the impact of LED can be sizable but only as a small 
perturbation on top of the standard oscillation unless we 
consider LED parameters not allowed by the disappear- 
ance modes. While LED can be potentially harmful in 
the appearance modes for the determination of the mass 
hierarchy and/or CP phase delta, we believe that the ap- 
pearance mode is not important (due to much smaller 
statistics than the disappearance ones) in constraining 
LED. 



III. CURRENT EXPERIMENTAL LIMITS 

Here we discuss the limits on the size of LED one can 
obtain from the former and current neutrino oscillation 
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experiments CHOOZ, KamLAND and MINOS. We could 
have considered other terrestrial experiments in our anal- 
ysis, but we have restricted ourselves to these three. Re- 
garding the long baseline experiments, KamLAND and 
MINOS are currently the best ones in terms of statistic 
and systematics. While the inclusion of other short base- 
line experiments could, in principle, improve our results, 
we have verified that this improvement is not very signif- 
icant, since a large fraction of the uncertainties of these 
experiments are correlated. 

We do not consider solar and atmospheric neutrino 
data for simplicity, and also because we do not expect 
significant improvement in constraining LED by adding 
these data (see Sec, |Y]). 



Reactor — > u e Experiments: CHOOZ and 
KamLAND 



The CHOOZ experiment is a former long (for reactor) 
baseline reactor neutrino oscillation experiment. Its goal 
was to probe the atmospheric oscillation parameters in 
order to shed light on the atmospheric anomaly [T5] . To 
achieve that aim, the experiment detected v e produced by 
the French CHOOZ nuclear power plant via the inverse 
/3-decay reaction D e + p — > e + + n. The 9 e energy E is 
estimated from the observed prompt energy E p of e + and 



nucleon mass difference M n 



M p as E 



E p + (M n 



M p ) + 0(E Pe /M n ), where the last term corresponds to 
the neutron recoil. Hence, the reaction has a 1.8 MeV 
threshold. 

The KamLAND (Kamioka Liquid scintillator Anti- 
Neutrino Detector) is a reactor neutrino oscillation exper- 
iment that operates in the site of the former Kamiokande 
experiment in Japan. Since 2003 KamLAND has ob- 
served T> e disappearance [11] compatible with the stan- 
dard neutrino oscillation scenario, giving strong support 
to the MSW LMA solution to the solar neutrino problem 
reported by the solar neutrino experiments [5j . The Kam- 
LAND detector observes v e produced by the surrounding 
nuclear power reactors via the same inverse /3-decay re- 
action described above. 

We have analyzed the last result by CHOOZ [19] and 
the most recent KamLAND data [H]. In fitting CHOOZ 
(KamLAND) data we have used the results of the new 
flux calculation for reactor neutrinos [27] [28] and varied 
Awijj and #13 (Am^ and #12) freely. We note, however, 
that the change of the reactor neutrino flux to the new 
one has very little impact on our results in obtaining LED 
bounds. For both experiments, we considered priors on 
all other standard oscillation parameters as explained in 
Appendix [B] except when comparing our standard Kam- 
LAND fit to [12, where we took 13 = 0. 

In Fig. [2] we show the region in the sin 2 26*13 — jAm^l 
plane allowed by CHOOZ data at 90% C. L. for the stan- 
dard oscillation case with a — 0. While we used the new 
reactor neutrino flux [27, 28 to study the impact of LED 
throughout this paper, in order to compare the results 



of our analysis with the original results by the CHOOZ 
group [19] (indicated by the solid red curve) we show the 
result obtained by using old flux (dashed blue curve) in 
addition to the one with the new flux (solid blue curve). 
We note that our simulation using the old reactor fluxes 
agrees reasonably well with that of CHOOZ [19 . 

We verified that the inclusion of LED does not essen- 
tially change this region. This can be understood if we 
remember that the main effect of LED is to induce os- 
cillations to the sterile KK modes. Since CHOOZ basi- 
cally does not see any significant deviation of the average 
v e — » v e probability from unity and the inclusion of LED 
can only lower this probability, LED cannot enlarge the 
CHOOZ allowed region. 

In Fig. [3] we show the regions in the tan 2 6*12 — Am^ 
plane allowed by the KamLAND data at 95%, 99% and 
99.73% C. L. for the standard oscillation case with a = 
(indicated by the dotted, dashed and solid curves) su- 
perimposed on the case fitted with LED (shaded colored 
regions). When we fit with LED we have also varied 
freely a, mo and the mass hierarchy. 

We see that our simulation agrees reasonably with the 
standard result of Fig. 2 of Ref. 12], our best fit point 
corresponds to Am^ = 7.84 x 10~ 5 eV 2 and tan 2 6*12 = 
0.46 with x r 2 „i„/dof = 17.3/15 = 1.15. With LED the 
allowed region gets considerably larger, our best fit point 
here corresponds to Ato^ = 8.28 x 10 -5 eV 2 , tan 2 612 = 
0.38 and a = 0.52 /im for NH with to = 3.56 x 10~ 2 eV. 
However, X 2 lin /dof = 16.8/13 = 1.29, so the inclusion of 
LED does not improve the fit. 
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CHOOZ 90% (hep-ex/0301017) 
CHOOZ 90% - old fluxes 
CHOOZ 90% - new fluxes 




sin 2 26» 13 

FIG. 2: Allowed regions in the sin 2 2^13 — lAm 2 ,^ plane ob- 
tained by fitting the CHOOZ data at 90% C. L. The result of 
our simulation using the old reactor fluxes is indicated by the 
dashed blue curve, which is to be compared with the result 
of the analysis A exclusion limit presented in [19] (red curve) . 
The solid blue curve represents the allowed region using the 
updated reactor fluxes from |28| . 
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KamLAND Standard KamLAND with LED 
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FIG. 3: Allowed regions in the tan 2 #12 — Am^i plane obtained 
by fitting the KamLAND data at 95%, 99% and 99.73% C. L. 
We compare the standard oscillation scheme (lines) with the 
LED oscillation scheme (colored regions). 
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FIG. 4: Excluded regions in the a — mo plane (mo is the 
lightest neutrino mass) by CHOOZ data at 90% and 99% C. 
L. for NH (blue curves) and IH (red curves). 



We have also investigated what region in the a — too 
plane can be excluded by CHOOZ and KamLAND data. 
This was calculated for NH and IH at 90% (99%) C. L., 
by imposing x 2 > Xmin + 4.61 (9.21), and is presented in 
Fig. |4j (CHOOZ) and Fig. [5] (KamLAND). As expected 
the v e — > v e channel gives a much more stringent limit 
on LED for the IH case (see FigJT}. We see that CHOOZ 
limits are stronger than KamLAND limits. For some 
numerical limits see Table HI 

We can understand qualitatively the shape of our ex- 
clusion curves in Figs. [4] and [5] as follows. If too > 0.05 
eV, neutrino masses are degenerate and in this case the 
limit has to be proportional to £ = \/2m a : i.e., if 
£ > Cmax the region is excluded; this explains the lin- 
ear behavior at the upper part of the plots of Figs. [4] and 
[5| If, however, too << 0.05 eV, LED will be constrained 
by £2,3 (NH) or £ 1)2 (IH) so the limit will not depend on 
TOq; this explains lower part of the plots. 
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FIG. 5: Same as Fig. [4] but excluded by KamLAND data. 



B. Accelerator — > Experiment: MINOS 

MINOS (Main Injector Neutrino Oscillation Search) is 
a neutrino oscillation experiment at Fermilab that has 
been running since the 2006 accelerator-beam disap- 
pearance [3 H] supporting the results from K2K [7] and 
the atmospheric neutrino experiments [6]. MINOS has a 
magnetized near detector with 29 t fiducial mass at 1.04 
km from the production target and a magnetized far de- 
tector with a fiducial mass of 4 kt at 735 km. Recently 
MINOS has also reported the observation of accelerator- 
beam disappearance |10j . which we will not consider 
in this work due to low statistics. In MINOS are 



identified by charged current interactions and the sign 
of the associated muon produced which is determined by 
the muon curvature under the detectors magnetic fields. 
The main background is due to neutral current events. 

We have analyzed the most recent MINOS data in the 
Vy, v n mode [TU] • In Fig. |6]we show the allowed regions 
in the sin 2 26> 23 - lAmfJ plane at 68% and 90% C. L. In 
the upper panel we have the pure standard oscillation (no 
large extra dimension allowed, a = 0) and in the lower 
panel we have allowed for LED. In fitting the data we 
have varied IAto^I and sin 2 2# 2 3 freely, and considered 
priors on all other standard oscillation parameters (see 
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FIG. 6: Allowed region for the standard oscillation parame- 
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|Am3i plane from MINOS 
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data. In the upper panel we assumed no LED while in the 
lower panel we allowed for LED in the fit. 



Appendix [B] for further details) . When we fit with LED 
we have also varied freely a, tuq and the mass hierarchy. 
Our best fit point in the standard oscillation fit is 
2.39 x 1(T 3 eV 2 and sin 2 26> 23 = 1 with 
/dof = 12.3/12 = 1.02. With LED the allowed re- 
gion gets enlarged, however the best fit point remains 
the same with a — 0, hence any value of mo is allowed. 
The Xmin/ dof = 12-3/10 = 1.23, so the inclusion of LED 
worsens the fit to data. 

We have investigated what region in the a — mo plane 
can be excluded by MINOS — > data. In Fig. [7 we 
present the excluded region calculated for NH and IH at 
90% and 99% C. L. As expected the — > channel 
is equally sensitive to NH and IH (see Figfl]). For some 
numerical limits, see Table [T] 



C. CHOOZ, KamLAND and MINOS Combined 



We have analyzed MINOS Vn — > Vu together with 
CHOOZ and KamLAND data by minimizing their added 
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FIG. 7: Same as Fig. [| but excluded by MINOS 
data. 



up x 2 functions letting all parameters vary freely. The 
excluded region for LED given by the combined fit is 
shown in Fig. [8] We see that the combined fit improves 
the limits derived until here, except for NH when m — > 
where the limit is basically that given by MINOS. For 
some numerical limits, see Table [T] 




FIG. 8: Same as Fig.|4]but excluded by CHOOZ, KamLAND 
and MINOS combined data. 



IV. FUTURE TERRESTRIAL NEUTRINO 
OSCILLATION EXPERIMENTS 

Here we discuss the possibility of improving the current 
limits on LED by the future neutrino oscillation experi- 
ments Double CHOOZ, NOi/A and T2K. 
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A. Reactor — > v e Experiment: Double CHOOZ 

The Double CHOOZ experiment [50], is a reactor neu- 
trino oscillation experiment that is being built in France 
which aims to explore the range 0.03 < sin 2 2#i3 < 0.2. 
There will be two identical 8.3 t liquid scintillator de- 
tectors, one at 400 m and the other at 1.05 km from 
the nuclear cores. The expected luminosity is 400 t GW 
y. We will consider 3 years of data taking in our calcu- 
lations. In fitting the data we have varied lAmfj-J and 
sin 2 2(9 13 freely, and considered priors on all other stan- 
dard parameters (See Appendix |b|) . 

In Fig. [9] we show our expected sensitivity for sin 2 2# 13 
as a function of (Amf^ for Double CHOOZ after 3 years 
for the standard oscillation analysis. We have verified 
that allowing for LED in the fit does not change this 
sensitivity curve as long as a < 0.3 /xm. So LED cannot 
simulate a nonzero #13. 



Double CHOOZ 




10" 2 10 _I 
sin 2 26 13 



FIG. 9: Sensitivity to sin 2 26 13 predicted for Double CHOOZ 
after 3 years, without assuming LED (standard oscillation). 
Here we have assumed as input: sin 2 2^13 = 0, lAm 2 .^ = 
2.46 x 1(T 3 eV 2 , and a = 0. 



We also have estimated the improvement that this ex- 
periment can provide on the limits given by CHOOZ and 
KamLAND. In Fig. [TO] we plot the potential exclusion 
region on the a — too plane As in the case of CHOOZ and 
KamLAND (see Figs. [4] and [5|, we obtained the better 
sensitivity for the IH case. For some numerical limits, 
see Table [H 
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FIG. 10: Sensitivity to LED predicted for Double CHOOZ 
after 3 years of data taking. 



B. Accelerator i/ M — > Experiments: T2K and 
NO^A 

T2K (Tokai to Kamioka) .21] is an experiment cur- 
rently running in Japan using a 0.75 MW beam 
from the J-PARC facility aimed at the 22.5 kt wa- 
ter Cherenkov detector Super-Kamiokande with a 295 
km baseline. T2K in its first phase will take data in 
Vy, -> mode. 

NOj/A (NuMI Off- Axis v e Appearance) [23], is an ex- 
periment that is currently being built in Fermilab and it 
will observe — > v^, e and £y, — > £^, e - The experiment 
will consist of a 222 t totally active scintillator detector 
(TASD) near detector and a 25 kt TASD far detector at 
810 km and 1.12 MW of beam power. 

We have simulated these experiments according to Ap- 
pendix [51 considering 5 and 3 years of — > data for 
T2K and NO^A, respectively. In fitting the data we have 
varied lAmjjjJ and sin 2 2# 2 3 freely, and considered priors 
on all other standard parameters. 

In Fig. [TT] we show the potential excluded region by 
T2K (5 yr) and NO^A (3 yr) in the a - m plane The 
limits are basically the same for those two experiments 
and they do not depend on the mass hierarchy. For some 
numerical limits, see Table [I] 

We see that neither of these experiments can really im- 
prove MINOS limits. The reason for that is the fact that 
LED induces oscillations into KK modes which are more 
sizable at higher energies away from the oscillation mini- 
mum (see Fig{T]) as the probability is larger in this region. 
T2K and NO^A are narrow (off-axis) beam experiments 
designed to measure precisely mixing parameters from 
the behaviors of oscillation probabilities around the first 
oscillation minimum, which means they are not very sen- 
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FIG. 11: Sensitivity to LED predicted for T2K (top panel) 
and NO^A (bottom panel) after 5 and 3 years of data, re- 
spectively. 



sitive away from it. 



V. DISCUSSION AND CONCLUSIONS 

We have investigated the effect of LED in neutrino os- 
cillation experiments assuming that singlet SM fcrmion 
fields can propagate in the bulk of a d-dimensional space- 
time and couple to the SM neutrino fields that lie in the 
brane through Yukawa couplings with the Higgs. We 
have shown that terrestrial neutrino oscillation experi- 
ments can provide submicrometer limits on the largest 
extra dimension a. 

For hierarchical neutrinos with m — > 0, CHOOZ, 
KamLAND and MINOS together constrain a < 
0.75(0.98) ^m for NH and a < 0.49(0.57) ^m at 90 (99)% 
C. L. for IH. For degenerate neutrinos with too = 0.2 eV 
their combined data constrain a < 0.10(0.12) fim at 90 
(99)% C. L. 





Limit on a ((im) at 90% (99%) C. L. 


Experiment 


NH, m -> 


IH, mo — > 


m = 0.2 eV 


CHOOZ 




0.54(0.61) 


0.13(0.14) 


KamLAND 




0.79(0.91) 


0.19(0.22) 


MINOS 


0.73(0.97) 


0.73(0.97) 


0.12(0.16) 


Combined 


0.75(0.98) 


0.49(0.57) 


0.10(0.12) 


Double CHOOZ 




0.38(0.46) 


0.09(0.11) 


T2K 


0.76(0.89) 


0.76(0.89) 


0.13(0.16) 


NOi^A 


0.80(0.92) 


0.80(0.92) 


0.14(0.17) 



TABLE I: Limits on the size a of the extra dimension for 
both hierarchies and degenerate neutrinos. See text for more 
details. 



We have also found that the future Double CHOOZ ex- 
periment will be able to improve these limits by roughly 
20% for the IH and 10% for the degenerate case. How- 
ever, T2K and NOi/A due to their narrow beam cannot 
surpass MINOS limits. 

Let us discuss briefly what we can expect from solar 
and atmospheric neutrinos. For solar neutrinos, if the 
ratio 1/a is much larger than yj Am^, the matter effect 
is not important as long as the impact of LED on the 
standard oscillation is concerned, and the effect of LED 
is to induce vacuum like oscillations from active to ster- 
ile states, simply reducing the overall v e (or all active 
v) survival probability. Since the inverse of the bound 
we obtained fr om MI NOS and KamLAND on a is much 
larger than \J Am^, we expect only a small reduction 
of solar v e due to LED which would not spoil the good- 
ness of fit of solar neutrinos by the standard MSW effect. 
In fact in order to induce strong distortion of the solar 
neutrino spectra, the size of a should be in the range of 
~ (60 — 100) fim [TS], much larger than the bound we 
obtained. Therefore, we expect that addition of the solar 
neutrino data to our analysis would not improve much, 
if at all, the bound on the size of the LED we obtained 
in this paper. 

For atmospheric neutrinos, we have checked that for 
given values of the size of the LED (a) and the light- 
est neutrino mass (too) and mass hierarchy, the magni- 
tude of the impact of LED on the — > (y^ — > v^) 
and v e — > v e (v e — > D e ) survival probabilities are similar 
to what we see in Fig. [T] for the relevant range of L/E 
from 1 — 10 4 km/GeV. This was done including earth 
matter effects making use of the formalism presented in 
Appendix [A] . We have verified that, as long as we con- 
sider parameters excluded by MINOS and/or KamLAND 
(shown in Figs. 4-6), LED does not make the oscillation 
probability deviate strongly from the standard oscillation 
scheme for atmospheric neutrinos. Therefore, by adding 
atmospheric neutrino data to our analysis, we do not ex- 
pect significant improvement on the bounds obtained on 
LED in our paper. 

Let us try to make some comparison of our bounds 
with the ones from the LHC. In LED models the con- 
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nection between the fundamental scale of gravity, Mjj, 
the number of extra dimensions d, the size of the com- 
pactification radius a and the Planck mass Mp is given 
by Mp +2 = Mp/(8Tra d ) 29J where it was assumed that, 
for simplicity, the size of the all LED radii a is equal 
for d > 2. So for d = 1 (d = 2) our limits on a imply 
M D > 10 6 TeV {M D > 22 TeV). On the other hand, 
the presence of LED also predicts gravition-emission and 
graviton exchange processes at colliders and according to 
Ref. [30] ATLAS and CMS at the LHC, after 36 pb" 1 , 
can exclude M D < 3 — 4 TeV, for d = 1 and 2, depending 
on the ratio A/Mc, A being the cutoff scale. Therefore, 
despite that the bounds we obtained in this work are 
model dependent, so far, they are stronger than the ones 
that come from collider physics. 

Recently, new flux calculations for reactor neutrinos 
became available [27] . We took them into account in our 
analysis of CHOOZ, KamLAND and Double CHOOZ but 
the impact of the change of the flux on our results is very 
small. Nevertheless, with this new flux calculation, older 
reactor neutrino oscillation experiments exhibit the so 
called reactor antineutrino anomaly recently reported in 
Ref. [25]. We note that, although the inclusion in our 
analysis of these older reactor neutrino oscillation exper- 
iments would not improve essentially the limits obtained 
in this work, they could favor some range of the LED pa- 
rameters currently allowed (obtained in this work), and 
therefore, deserve further study [31 . 

A final comment is in order. One cannot directly apply 
our limits to models such as the one discussed in Ref. [35] , 
where neutrino oscillations are modified by the presence 
of reconstructed nongravitational large extra dimensions. 
However, since the model studied here is the continuum 
limit of the former, we suspect that similar constraints 
could be derived in that case. 



with respect to the active flavors by defining the unitary 
transformations 
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so that Y^af) Uai m af} Rpj = ^ij Mi , with the lower case 
Roman indices i,j = 1,2,3. Throughout this paper 
Greek indices will run over the 3 active flavors, Roman 
lower case indices over the 3 SM families and upper case 
Roman indices over the KK modes. Explicitly 
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where £j = \f2rrii a. 

Let us define the following states 
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Appendix A: Solution of the Evolution Equation for 
a Constant Matter Potential 

While matter effects are not very important for this 
work, in this appendix, for the sake of completeness, we 
describe the solution of the evolution equation in the 
presence of constant matter potential in the context of 
large extra dimensions. See also [161 118j where a sim- 
ilar procedure was adopted. We first diagonalize rn®a 



^(^M%f...) T , i = l,2,3 (A4) 

so that 
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This allows us to write the neutrino evolution equation 
in matter as 
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(A7) 



where £7 is the neutrino energy and we have defined 



V a \ I S ea Vcc + V NC 
) \ 



(A8) 



with the matter potentials Vcc = \/2G F n e and y NC = v< p ) ^ L ) in Ec l- fs} 



electron (neutron) number density in the medium and 
NC stands for neutral current. 

Here we describe how to obtain an analytic expres- 
sion for the eigenvalues and the amplitudes 
needed to calculate the transition amplitudes A{va^ — > 



— ^Gpn n . Gf is the Fermi constant, n e (n n ) is the 



If we multiply Eq. ( A2 ) by its conjugate we get 
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where 



Vi = (N+l/2)i 
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i ' 



(A10) 



= 1,2,3 the generation indices. 

Defining Vy = 2£a 2 ^ U* m U aj V a we can reorga 



Ui = (6 2& • • ■ 2V&) and if = diag (l 4 9 ... A 2 ) with nizc Eq. flA7| 
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To diagonalize % we have to find the eigenvalues A, that solve det (2Ea 2 H — A 2 /) = 0. One can show, by 
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using the Gauss algorithm for determinant calculation To find the eigenvectors wf , corresponding to the 
that this is equivalent to calculate eigenvalues \\ N} we have to solve 



det (T) = 0, 
where T is a 3 by 3 matrix with elements 



(A12) 



N _ \{N)2 N 



H wf = A 



(A14) 



T- ■ — 



-A 2 + ^ cot (ttA) 



Sij + Vij, = 1,2,3). 



(A13) 



where we denote an element of wf by {wf )| = W i 



,N\M - t ff (™) 



Explicitly in terms of these elements Eq. (A14| can be 
written as: 



K 



(A15) 



A=l 



1=1 



and 



We can obtain an equation for W^f from Eqs. (|A15|), 



(^-(Af'f)^'^. (A16) 



and Eq. (A16) and Eq. (A10) in the limit N -> oo: 
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(A17) 



* < 0) ( cot (.A^) - (Af >) 2 ) + E ^ 
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(A18) 



So that for each eigenvalue obtained by solving to obtain W^f 10 ^ . We also need to impose the normaliza- 



Eq. (A12| one has to solve 

3 



tion of the eigenvector with 



Ewf 0) = o, 



(A19) 
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i=i i 
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= 1. 



(A20) 



As a technical note: in practice it is a very good ap- 
proximation to consider only the first five KK modes 
in the numerical calculation. We have verified that the 
inclusion of higher modes do not cause any significant 
change in our results. 

In vacuum, T tJ = T { and wjf 0) = V^ (JV0) as the KK 



modes connected to different generations decouple. In 
this case, if a -1 << m i? as show in Ref. [15] , we have 



(wr ] j 



1 + 0(g) iV = 
+ N = 1,2, 3. 



(A21) 
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Appendix B: Simulation Details 

In this section we gather all information used to sim- 
ulate the experiments. We implemented all experiments 
using a modified version of GLoBES 33J . To model the 
energy resolution, we used a following Gaussian smearing 
function, 



R (E, E') = 



1 





(E-E 1 ) 2 

' — ^ — 



(Bl) 



where <te was defined according to each experiment (see 
below). 

Let us name the % 2 function without any uncertainty 
and previous knowledge of oscillation parameters as Xo- 
To account for previous knowledge on some set of oscil- 
lation parameters we use Gaussian priors. Consider that 
these parameters pi have mean values pi and mean devi- 
ations Upi. Then, the Gaussian priors are added to the 



X 2 as 



2 2 1 
X = Xo + 



E 



{pi -PiY 



(B2) 



pi 



To deal with an experimental uncertainty (in flux, fidu- 
cial mass, etc), we modify Xo ~~ > Xo by adding a new pa- 
rameter x and add a penalty term x 2 /er 2 . To exemplify 
that, let us assume an uncertainty ctnc m the neutral cur- 
rent events normalization. If A 4 NC is the number of neu- 
tral current events simulated at the «-th bin, then in the 
Xq function we could replace N^ c — > (1 + £nc)^ NC an d 
add the penalty term £nc/ nc to the resulting x 2 func- 
tion. In summary, taking into account previous knowl- 
edge in the oscillation parameters and experimental un- 
certainties, the resulting x 2 has the form 



calculated number of events in the i-th energy bin which 
depend on the standard oscillation parameters and, in 
the case of LED, also on mo, a and the neutrino mass 
hierarchy. Given the complexity of the inclusion of mat- 
ter effects in the LED framework (see Appendix |A| , our 
simulations were done using the vacuum oscillation prob- 
abilities. For KamLAND and MINOS this is acceptable 
because the matter effects play a small role on the sur- 
vival channels. For CHOOZ and Double CHOOZ, since 
the baseline is short, the matter effects are negligible. 
Finally for T2K and NO^A, as long as we arc fitting 
simulated data, the matter effects are important only in 
the appearance channels, which are not used. 



1. CHOOZ 

In order to reproduce Fig. 55 of Ref. [15] , we con- 
sidered a detector located at 1.05 km from the nuclear 
cores. The predicted antineutrino spectrum was based on 
the newest fluxes calculation available [37] and the over- 
all normalization was chosen so that the ratio between 
the observed and theoretical unoscillated total number 
of events would match the value given by |28j . which is 
R™p = 0.961. 

Our analysis was based on rates information only, so 
we minimized a % 2 function composed by 



Xo 



/i? cx P - R 



thco 



(B4) 



with respect to all parameters considered free in the fit. 
Here R theo is the ratio between the observed and theoret- 
ical oscillated total number of events and a = 4.2% takes 
into account the statistical and systematical uncertainty. 



xl 



(pi-PiY 



E 

3 



j 

r 2 . 

X] 



(B3) 



For a detailed explanation about these techniques, see 
the GLoBES manual [55] , 

Generically, for the data fits we have varied some of 
both standard and LED oscillation parameters. When 
we considered Gaussian priors for a standard parameter, 
we based the previous knowledge on [33], using, at la, 
Amlj = 7.59±0.20xl0~ 5 eV 2 , Amli = 2.46±0.12x 10" 3 
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eV 2 , 6>i2 = 34.4° ± 1°, 6*23 = 42.8° ± 4.7° and we used a 
conservative limit for (9 13 , sin 2 20 13 < 0.09. We did not 
impose any prior on 5c p- For all fits with LED we varied 
freely a, mo and the mass hierarchy. 

It is useful to define the following quantities before giv- 
ing the details of each experiment. For KamLAND and 
MINOS, N° xp are the experimental data points taken 
from Fig.l of Ref. [12] and Fig. 1 of Ref. [10], respectively, 
while for the future experiments JV™ P are the simulated 
data points calculated assuming fixed values for oscilla- 
tion parameters. Moreover, iV* heo are the theoretically 



2. KamLAND 

We follow our previous papers [35] in calculating the 
number of events expected from reactors for a total ex- 
posure of 2881 t yr. However, to calculate the unoscil- 
lated v e spectrum we have updated the averaged ratios 
of the fission yields of the four isotopes that significantly 
contribute to the flux as: 235 U: 238 U: 239 Pu: 241 Pu = 
0.570: 0.078: 0.295: 0.057, in accordance with Ref. [12] . 
The energy resolution was modeled as a Gaussian with 
a E = 0.064 y/E/MeV- 0.8. 

We have determined the experimentally allowed re- 
gions minimizing the x 2 function composed by 
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x^ = E 

8=1 



(iVf p - Nl heo ) 2 



cxp2 ' 



(B5) 



with respect to all parameters considered free in the fit. 



Here a s 



4.3%. The experimental data points were 
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taken from Ref. [T3] in the energy window from 1.7 to 
8.925 MeV (bin width of 0.425 MeV). All uncertainty is 
included in cr sys and we used the efficiency given in [T2"] . 



3. MINOS 

MINOS simulation was performed in accordance with 
[55] . using the NuMI neutrino beam given by [37], the 
neutrino-nucleon cross section from [351 ISH] • The analysis 
was performed with neutrinos in 250 MeV bins from 1 to 
5 GeV. We assumed uncertainties in the signal and back- 
ground that were taken to be er s = 4% and o^c = 3%, 
respectively. The detecting efficiency was taken from 
Ref. [lOj and the energy resolution was modeled as a 
Gaussian with a E = 0.16 £7GeV+0.07 y/E/GcV to best 
reproduce the MINOS allowed region for the standard os- 
cillation parameters. 

We have determined the experimentally allowed re- 
gions minimizing the x 2 function composed by 
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Xo 



\iV* heo 



(B6) 



with respect to all parameters considered free in the fit. 



Double CHOOZ 



by a Gaussian with <je = 80 MeV for signal reconstruc- 
tion and 2% uncertainty in the flux and background. For 
more details see [4Tj . 

To estimate T2K sensitivity we minimize the \ 2 func- 
tion composed by 
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2 _ (JV^P-iVjheo)S 



i=l 



(B8) 



with respect to all parameters considered free in the fit. 



6. NOi/A 

The experimental setup considered was based on 
[231 HQ], being a 25 kt TASD far detector at 810 km, 
1.12 MW of beam power, 3 years of data taking in 
the Vfj, — > mode, 20 energy bins from 1 GeV to 
3.5 GeV and energy resolution modeled by a Gaussian 
with oe = 0.05 ^E/GeV for signal reconstruction and 
&e = 0.10 WEjGeV for neutral current reconstruction. 
We assumed uncertainties in the signal and background 
normalization using a slightly different method as dis- 
cussed above. We used method "C" of GLoBES man- 
ual [33] with a and b parameters (5%:2.5%) for both sig- 
nal and background. 

To estimate NOjaA sensitivity we minimize the x 2 func- 
tion composed by 



Basing Double CHOOZ simulation on j20||40j, we used 
two identical 8.3 t liquid scintillator detectors, one at 400 
m and the at 1.05 km from the nuclear cores. The ex- 
pected luminosity is 400 t GW y. We considered 3 years 
of data taking assuming 62 energy bins from 1.8 to 8 MeV 
with the energy resolution modeled by a Gaussian with 
a E = 0.12 yjE/MeV — 0.8. The uncertainties taken into 
account for both cores and detectors were: isotopic abun- 
dance (2%), core power (2%), flux normalization (0.6%), 
overall flux normalization (2.5%) and energy scale for 
each core (0.5%). 

To estimate Double CHOOZ sensitivity we minimize 
the x 2 function composed by 



2(1 



2 _ (Nr _ N t^ 
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(B9) 



with respect to all parameters considered free in the fit. 
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i=l d 
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rCXp2 • 
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(B7) 



where a sys — 1%, with respect to all parameters consid- 
ered free in the fit. 



5. T2K 



We base T2K simulation on Ref. [H] where we have 
considered a beam power of 0.75 MW, a 22.5 kt water 
Cherenkov detector at 295 km from the neutrino source, 5 
years of data taking in the v„ — ¥ z/„ mode, 36 energy bins 
from 0.2 GeV to 2.0 GeV and energy resolution modeled 
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